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Abstract 

The electronic structure of the low-temperature (LT) monoclinic magnetite, Fe304, is investi- 
gated using the local spin density approximation (LSDA) and the hSDA+U method. The self- 
consistent charge ordered hSDA+U solution has a pronounced [001] charge density wave character. 
In addition, a minor [00|] modulation in the phase of the charge order (CO) also occurs. While 
the existence of CO is evidenced by the large difference between the occupancies of the minority 
spin t2g states of "2-|-" and "3-|-" Fe^ cations, the total 3d charge disproportion is small, in accord 
with the valence-bond-sum analysis of structural data. Weak Fe orbital moments of ~0.07/iB ^.re 
obtained from relativistic calculations for the CO phase which is in good agreement with recent 
x-ray magnetic circular dichroism measurements. Optical, magneto-optical Kerr effect, and O K- 
edge x-ray absorption spectra calculated for the charge ordered LSDA+U solution are compared 
to corresponding LSDA spectra and to available experimental data. The reasonably good agree- 
ment between the theoretical and experimental spectra supports the relevance of the CO solution 
obtained for the monoclinic LT phase. The results of calculations of effective exchange coupling 
constants between Fe spin magnetic moments are also presented. 

PACS numbers: 7L20.-b, 7L28.-Kd, 7L30.+h 
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I. INTRODUCTION 



The problem of a theoretical description of metal-insulator transitions has a challenging 
history of almost 70 years. It was first addressed by Verwey, de Boer, and Peierls in the 
late 1930's; they pointed out the extremely important role of electron-electron correlations 
in a partially filled (i-electron band in transition metal oxides,-"^ e.g., nickel oxide (NiO) 
and magnetite (Fe304). In both systems the metal-insulator transition occurs, violating the 
Bloch- Wilson band-insulator concept, the only one known at that time.-'^'^'^ These earlier 
observations launched the long and continuing history of the field of strongly correlated 
electrons. In the past 70 years, much progress has been achieved from both theoretical 
and experimental sides in understanding strongly correlated electrons and metal-insulator 
transitions.- However, the charge ordering proposed by Verwey behind the metal-insulator 
transit ion^i^iii^ in Fe304 remains at the focus of active debate.— 

Magnetite is a permanent natural magnet. Its magnetic properties have fascinated 
mankind for several thousand years already.— Comprehensive investigation of magnetite 
started in 1913 when the first indication of the first-order transition was obtained by a sus- 
ceptibility measurement on a synthetic polycrystal.— An anomalous transition peak at 
~130 K was found in addition to the already known paramagnetic-to-ferrimagnetic order 
transition on cooling through ~850 K. Electric resistivity measurements initially performed 
by OkamuraW- and Verwey- showed that this anomalous behavior of the susceptibility co- 
incides with a sharp two orders of magnitude increase of resistivity on cooling below ~120 
K.—^ The formal valency of Fe^ cations on the octahedral B sublattice of the inverted spinel 
AB2O4 structure is noninteger (2.5+). Verwey suggested that this transition is caused by 
the ordering of Fes cations, with a simple charge arrangement of (001) planes alternately 
occupied by 2+ and 3+ Fe^ cations.-*^*^ This particular (Verwey) charge order (CO) obeys 
the so-called Anderson criterion^^ of minimal electrostatic repulsion leading to a short range 
CO pattern with each tetrahedron formed by B sites being occupied by equal numbers of 
2+ and 3+ cations. 

However, further experiments disproved the orthorhombic Verwey CO model. These ex- 
periments have clearly established the rhombohedral distortions of the cubic unit cell first 
detected by Rooksby et al. from x-ray powder diffraction.—!^ Furthermore, observations of 
superstructure reflections revealed half-integer satellite reflections, indexed as (/i. A;, / + 1) on 
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the cubic unit cell, which points to a doubling of the unit cell along the c axis and shows the 
symmetry to be monoclinic— The observation of monoclinic lattice symmetry was also 
confirmed by a single- crystal x-ray study,^'^ whereas the observation of a magnetoelectric ef- 
fect indicated even lower PI symmetry in the low-temperature phase. Although, clear evi- 
dence of the monoclinic lattice symmetry below Ty was obtained, small atomic displacements 
have not been fully resolved so far. The absence of a definitive, experimentally determined 
structure gives rise to many theoretical models proposed for the low-temperature (LT) phase 
of magnetite.— In particular, purely electronio2^i2ii2^i2£ and electron-phonon^SiSii^^ models for 
CO, as well as a bond dimerized ground state without charge separation,— have been pro- 
posed. 

At room temperature Fe304 is a poor metal with an electronic resistivity of 4 mQ cm 
which is considerably higher than the resistivity of simple metals (for instance, the resistivity 
of Cu is 1.7 fiQ cm). The microscopic origin of the electronic state above the Verwey 
transition is still unclear. Photoemission studies performed by Chainani et a/.^^ and Park et 
alr^ gave controversial results indicating metallic and semiconducting behavior above Ty, 
respectively. In accordance with Park et al.,— qualitatively similar conclusions have been 
derived by Schrupp et al.^^ and Pimenov et al?^ basing on a soft x-ray photoemission study 
and complex conductivity measurements at terahertz/infrared frequencies, respectively. 

Magnetite is a ferrimagnet with anomalous high Curie temperature of ~850 K. The A- 
site magnetic moments are aligned antiparallel to the i?-site moments. This remarkable 
situation of antiferromagnetic coupling between Fe^^ and Fe^^"'' cations, which have 
and d^'^ electron configurations in the high-spin state, respectively, should lead to a half- 
metallic state with integer magnetic moment of 4//^ per formula unit. This is in qualitatively 
good agreement with magnetization measurements, which result in a ~4.1/iB magnetic 
moment per unit cell.— In contrast, spin resolved photoemission measurements show ~20% 
suppression of the pure negative spin polarization from the half-metallic behavior obtained 
for Fe304 (111) films. Scanning tunneling microscopy and spin-polarized photoelectron 
spectroscopy give only 55% of negative spin polarization at the Fermi level.'^^ To explain 
this phenomenon two possible scenarios have been recently proposed. According to the 
first, Fe304 is still a half metal in the bulk, whereas surface stresses lead to unusual lattice 
distortions which reduce spin polarization at the Fermi level to ~40%.— An alternative 
explanation is the formation of nonquasiparticle or incoherent states in the gap near the 
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Fermi level due to correlation effects.— 

Recent bond- valence-sum analysis^ of high-resolution neutron and x-ray powder diffrac- 
tion data results in a small charge disproportion of only 0.2 e between Fe^ cations with the 2+ 
and 3+ formal valency.— This interpretation has been the subject of much controversy.— 
However, the smallness of the charge order parameter was reproduced in an electronic struc- 
ture study of the refined low-temperature crystal structure using the load spin density 
approximation (LSDA)+f/ method.— In particular, a more complicated charge-ordering 
pattern inconsistent with the Verwey CO model was obtained. In addition to that, the 
t2g occupancy self-consistently obtained in the LSDA+[/ calculations is strongly modulated 
between the Fe^"*" and Fe^"^ cations, yielding a distinct orbital order with an order parame- 
ter that reaches 70% of the ideal value.^ Since no direct experimental confirmation of this 
charge and orbital order pattern is so far available, the interpretation of these results is still 
open to debate. However, this behavior seems to be universal and has recently been found 
in several other charge ordered mixed-valent systems.— i^i^i^i^ 

In order to further check the pertinence of the CO model obtained self-consistently in 



Ref. 47| we carried out a detailed theoretical study of exchange coupling constants, optical 
conductivity, magneto-optical (MO) Kerr effect, and x-ray absorption at the O K edge 
of low-temperature Fe304 and compared the results of the calculations to the available 
experimental data. 

The paper is organized as follows. In Sec. [TTl we discuss the low temperature crystal 
structure of Fe304. Section [Till presents computational details relevant to reproduce the 
calculation results. The results of electronic structure calculations obtained by the LSDA 
and LSDA+f/ methods for the low-temperature Fe304 are presented in Sees. [IV] and |Vl 
respectively. In Section |V] a charge and orbital order are also discussed and results of cal- 
culations of exchange coupling constants are presented. In Section |Vl] we present calculated 
optical, magneto-optical, and O K-edge x-ray absorption spectra and compare them to 
corresponding experimental spectra. Finally, the results are summarized in Section IVIIi 

II. CRYSTAL STRUCTURE AND CHARGE ORDER 

Above the Verwey transition magnetite crystallizes in the face-centered cubic (fee) inverse 
spinel crystal structure with space group Fd3m. The iron atoms occupy the interstitial 
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positions of the close-packed fee strueture formed by the oxygen atoms. According to the 
chemical formula of inverse spinel, AB20i, there are two different iron sublattices {A and B) 
distinguished by the point symmetry (T^ and D^d) as well as by the averaged Fe-0 distance 
(1.876 and 2.066 A for A and B iron sites, respectively) and the valence state of iron cations. 
In particular, the A sublattice is formed by Fe^^ cations tetrahedrally coordinated by four 
oxygen ions, whereas the so-called B sublattice consists of iron sites octahedrally coordinated 
by six oxygen ions. The octahedral -B-sites are occupied by an equal number of randomly 
distributed 2+ and 3-|- Fe cations, which results in an average valence value of 2.5+ per 
each Fes cation. The -B-sublattice is highly frustrated and can be considered as a diamond 
lattice of Fe^ cation tetrahedra, sharing corners with each other. 

One of the driving forces for formation of the charge-ordered state in mixed valent tran- 
sition metal oxides is Coulomb repulsion. In particular, another key to understanding the 
charge-ordered structure in Fe304 is provided by the Anderson criterion for the minimum 
of electrostatic repulsion energy.— According to Anderson each tetrahedron formed by octa- 
hedral B sites of the spinel structure is occupied by two 2+ and two 3+ Fe cations in order 
to minimize the intersite electrostatic repulsion. Thus, within the Cc super cell of the cu- 
bic Fe304, there are ten independent cationic arrangements, whereas one of them coincides 
with the Verwey model.^ However, the Verwey CO pattern has overall the lowest classical 
electron correlation energy. 

The low-temperature structure was shown to have a \/2a x ^/2a x 2a supercell with space 
group Cc from x-ray and neutron diffraction.— However, recent structural refinement (at 
90 K) was only possible in the centric monoclinic space group P2/c with ^ x ^ x 2a of the 
cubic spinel subcell and eight formula units in the primitive unit cell.— Since the refinement 
for the P2/c space group was found to be unstable, additional Pmca orthorhombic symmetry 
constraints were also applied. This is equivalent to averaging the true superstructure over 
the additional symmetry operators, i.e., each B site in the P2/c unit cell is averaged over four 
nonequivalent subsites in the large v^a x ^/2a x 2a Cc supercell. Note, however, that such 
an approximation is robust in the sense of smallness of any distortions from the P2/c subcell 
to the Cc monoclinic cell (according to Ref. |45| these are of ~0.01 A). Previous structure 
refinement below Ty obtained by lizumi et al. resulted in a x x 2a subcell of the 
Cc unit cell and imposed orthorhombic symmetry constraints on the atomic positions.— In 
particular, a refinement based on an approximation of the true crystal structure by a centric 
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space group Pmca or polar Pmc2i was proposed. But a charge-ordered arrangement has not 
been identified in this refinement, although large atomic displacements of Fe and O atoms 
were found. This is in strong qualitative contrast to the recent structure refinement proposed 
by Wright et al. where clear evidence of CO below the transition has been found.— 

According to the refinement the octahedral Fe^ sites are split into two groups with 
different values of the averaged Fe-0 bond distances; with B2 and -B3 sites being significantly 
smaller than Bl and BA [Bl-BA are crystallographically independent Fe^ sites according 
to the notation in Refs. lulandQ). A different averaged Fe-0 bond distance is a sensitive 
experimental indicator of the cation charge state. Quantitative analysis of the valence state 
of both Fe^ groups using the bond-valence-sum (BVS) method shows that the octahedral 
Fee sites fall into two clear groups with respect to the estimated value of valence. The 
result is a charge disproportion of 0.2 e between large {Bl and BA) and small {B2 and 53) 
sites (which has been referred as the class I CO model). Another possible class of CO arises 
from the symmetry-averaging orthorhombic constraint. There are 32 charge-ordered models 
which are refered to as class II CO because large (-B1 and 54) and small {B2 and 53) sites 
could be averaged over (3 Fe2+ + Fe^+) and (Fe^^ + 3 ¥e^+) subsites, respectively. The 
symmetry averaging results in decrease of the more pronounced charge separation of 0.4e in 
the full Cc superstructure (class II CO) down to 0.2 e in the P2/c subcell. The Anderson 
criterion is not satisfied by any of the class I or class II CO models. This is remarkable 
because the Anderson criterion has been widely used in many CO models.— However, 
class II, as was shown from electrostatic repulsion energy estimations, appears to be more 
plausible than the class I arrangement. 

Recently this interpretation of the refined crystal structure has been found to be contro- 
versial. The lack of atomic long-range CO and, as a result, an intermediate valence regime 
below the Verwey transition were proposed.— It is argued that the difference of the av- 
erage Fe-0 distances between compressed and expanded FeOe octahedra, which could be 
considered as a maximum limit of charge disproportionation, has the same order as the total 
sensitivity (including experimental errors) of the bond-valence-sum method. This remark- 
able controversy shows that the understanding of the system is far from satisfactory. 
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III. COMPUTATIONAL DETAILS 



In this paper we present results of band structure calculations which have been performed 
for the low-temperature crystal structure of Fe304 recently refined at 90 K.— The refined 
cell parameters used in calculations are a = 5.94437 A, b = 5.92471 A, c = 16.77512 A, 
and (3 = 90.236°. The monoclinic P2/c unit cell contains eight formula units. There are 
eight different types of iron atom sites: two tetrahedrally coordinated A sites and six B sites 
octahedrally coordinated by six O atoms. 

The electronic structure of magnetite was calculated self-consistently using the local- 
spin-density approximation and LSDA+f/ approacb^*^ with the linear muffin-tin orbitals 
(LMTO) method in the atomic-sphere approximation.— The radii of the muffin-tin spheres 
were taken as -RFe= 2.125 a.u. and -Ro= 2.0 a.u. Fifteen kinds of empty spheres were 
introduced to fill up the interatomic space. As was shown in Ref . |57|, weak spin-orbit coupling 
does not appreciably affect the band structure of Fe304. We neglect it for simplicity in Sees. 
HVl and |Vl Optical, magneto-optical (MO), and x-ray absorption spectra (XAS) for the 
P2/c model of the LT phase of Fe304 were calculated taking into account the spin-orbit 
coupling using the spin-polarized relativistic LMTO method.— In order to calculate O K- 
edge XAS spectra and the absorption parts of the optical conductivity in a wide energy 
range Fe 4/ and O 3d states were included into the LMTO basis set. The dispersive parts 
of the conductivity were calculated using the Kramers-Kronig relations. XAS spectra were 
calculated in the dipole approximation neglecting relaxation effects caused by the O Is 
core hole. A detailed description of the formalism for the calculations of MO and x-ray 
absorption spectra and corresponding matrix elements using the relativistic LMTO method 



can be found in Refs. 



57 



5S 



60 and references therein. 



IV. LSDA BAND STRUCTURE 

Figure [1] shows the total density of states (DOS) and band structure of Fe304 obtained 
from the LSDA calculations. Interestingly, the crystal structure distortion taken explicitly in 
the monoclinic phase does not strongly affect the electronic structure of Fe304. In particular, 
the calculation results agree well with previous band structure calculations for the cubic 
phase.— >^ Thus, the LSDA gives a uniform half-metallic ferrimagnetic solution with partially 
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FIG. 1: (Color online) Total density of states and band structure of the P2/c phase of Fe304 
self-consistently obtained using the LSDA method. The Fermi level is shown by dotted lines. The 
energy bands predominantly originating from the Fe^ 3d states are shown in light gray (green) 
color, whereas Fe^ t2g and eg bands are in dark gray (red) and gray (blue), respectively. The 
corresponding contributions to the total DOS are shown in the right panel. 

filled bands originating from the minority spin t2g orbitals of Fe^ cations. An energy gap of 
~1 eV opens in the majority spin channel between occupied Fe^ Cg and the bottom of the 
empty Fe^ states. The lower part of the valence band (below —3.5 eV) is mainly formed by 
O 2p states with a bonding hybridization with Fe 3d states, whereas bands near the Fermi 
level, between —3.5 and 2 eV, have a predominant contribution of Fe 3d states. 

It is worth recalling that Fe^ 3d states are split by the cubic component of a ligand field 
into a triplet t2g and a doublet Cg. Already in the Fd3m phase the local symmetry of Fe^ 
sites {D3d) is lower than cubic and the trigonal component of the ligand field splits states 
into a singlet aig and doublet e'g. In the monoclinic LT phase the symmetry is further lowered 
by the distortions and all the degeneracy is lifted. Nevertheless, the cubic component of the 
ligand field, which is determined by the relative strength of Fe 3d-0 2p hybridization of vr 
and a type, remains dominant, whereas the splitting within t2g and Cg subbands is smaller 
than the corresponding bandwidth. This allows one to label the corresponding states as t2g 
and Cg. Note, however, that the P2/c frame is rotated by ~ ^ with respect to the Fd3m one 
and the angular dependence of the t2g states is given by d^z ± dy^ and (i^.2_y2 combinations 
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FIG. 2: (Color online) Partial DOS obtained from the LSDA calculations for the P2/c phase of 
Fe304. The Fermi level is shown by dotted lines. 

of cubic harmonics. 

Figure [2] presents different contributions to the total density of states given by independent 
iron and oxygen sites. Crystallographically independent iron sites are labeled according to 



the notation of Ref. 
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Note, however, that the DOSs for two pairs of sites that were 
constrained to have equivalent coordination during the structural refinement - Fe^ia, Fe^ib 
and FeB2a; FeB2b ^ are very similar and we present the DOS averaged over the corresponding 
pairs, dropping a and b indices. 

The splitting of Fe^ states due to the octahedral component of the ligand field is roughly 
of 1 eV. This is considerably smaller than the exchange splitting between minority and 
majority spin states, which is ~3 eV, and consistent with the high-spin state of Fe cations. 
The absolute values of magnetic moments obtained by LSDA are 3.14 fis, 3.30 /i^, 3.44 /i^, 
3.27/iB, and 3.38 /i^ for Fe^ and FeBi-_B4 cations, respectively. The total ferromagnetic 
moment per Fe304 formula unit is 4.0 /ig. 

It should be noted that in contrast to experimental data the LSDA calculations result 
in a half-metallic solution without charge ordering. Apparently, the change of the LSDA 
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electronic structure, produced by the crystal structure distortion from the cubic to the mon- 
oclinic phase, is not sufficient to explain the charge ordering and metal-insulator transition 
in Fe304. The electron-electron correlations, mainly in the 3d shell of Fe cations, play a 
significant role. 

V. LSDA-ht/ RESULTS 
A. Band structure 

In order to account for the strong electronic correlations in the Fe 3d shell, at least 
on the static Hartree-Fock level, we calculated the electronic structure of the LT phase 
of Fe304 using the LSDA+t/ method. Previous studies proved this method to be quite 
successful in treating transition metal oxides with strong electron-electron correlations as 
well as systems with long-range order.— i^*^!^ The value of the U parameter for Fe cations 
estimated using different experimental and theoretical technics lies in the range of 4.5-6 
eV.— 1^1^ A reasonably good agreement of the calculated gap value of 0.18 eV with the 
experimental valued of 0.14 eV at 10 K was obtained using the U value of 5 eV. Note, 
however, that the charge and orbital order derived from the LSDA+f/ calculations does not 
depend on the exact U value when it is varied within the above mentioned limits. The value 
of the Hund's coupling J=l eV was estimated from constrained LDA calculations.— In the 
following all results presented in the paper were obtained using a U value of 5 eV. 

Figure [3] shows the LSDA+f/ band structure and the total DOS calculated self- 
consistently for the low-temperature structure of Fe304 using the Coulomb interaction pa- 
rameter U=5 eV and exchange coupling J=l eV. The corresponding partial Fe^ 3d DOS 
are shown in Fig. |H The LSDA+f/ calculations give results qualitatively distinct from those 
of the LSDA. An indirect energy gap of 0.18 eV opens in the minority spin channel between 
M and F symmetry points. One of the minority spin states of Fe^i and Fbba ions be- 
comes occupied while the FeB2 and Fe^s t2gj, states are pushed above the chemical potential. 
Although, as will be discussed below, the calculated disproportion of Fe^ 3d charges is sig- 
nificantly less than 1, in the following we use the notations Fe^"*" and Fe^"*" for FeBi,B4 and 
FeB2,BA cations, respectively, having in mind the difference of their t2g occupations. The 
top of the valence band is formed by the occupied t2gj, states of 51 and 54 Fe^"*" cations. 



10 



minority spin 




80 60 40 20 G X M Y G D 



FIG. 3: (Color online) Total DOS and band structure of the P2/c phase of Fe304 self-consistently 
obtained by the LSDA+C/ with U =5 eV and J=l eV. The Fermi level is denoted by the horizontal 
line and is taken as the zero of energy. A band gap of 0.18 eV opens between M and T symmetry 
points. The energy bands predominantly originating from the Fe^ t2g and eg states are shown in 
dark gray (red) and light gray (green) colors, respectively, whereas the gray (blue) color corresponds 
to Fe^"*" t2g bands. For the majority spin an energy gap of ~2 eV opens between Fe^"^ eg and Fe^ 
bands shown in dark gray (red) and light gray (green) color, respectively. The corresponding 
contributions to the total DOS are shown in the left panel. 

The bottom of the conduction band is formed predominantly by the empty t2gi states of 
B2 and 53 ¥e^+ cations. The remaining unoccupied t2gi states of Bl and BA Fe^"*" cations 
are pushed by the strong Coulomb repulsion to energies above 2.5 eV. Majority spin "id Fe^ 
states are shifted below O 2p states, which form a wide band in the energy interval between 
—7 and —2 eV. This is in strong contrast to the uniform half-metallic solution obtained by 
the LSDA. 

Bands corresponding to the Fe^"^ cations are fully occupied (empty) for minority (major- 
ity) spin states, respectively, and already in the LSDA do not participate in the formation 
of bands near the Fermi level. The LSDA+f/ method does not strongly affect these bands, 
which lie in the energy interval of —6 eV below and 1-2 eV above the Fermi level. 
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FIG. 4: (Color online) Partial DOS obtained from the LSDA+C/ calculations with U=5 eV and 
J=l eV for the P2/c phase of Fe304. The Fermi level shown by dotted lines. A charge-ordered 
insulating solution is obtained. Fe 3d minority states corresponding to SI and i?4 sites are occupied 
(Fe^"*") and located just below the Fermi level, whereas B2 and B3 are empty (Fe^^). The charge 
difference between 2+ and 3+ Fe^ cations is found to be 0.23 e. 



B. Charge ordering 

The obtained solution for CO of Fe^^ and Fe^"'' cations on the B sublattice is described by 
a dominant [001]c charge (and spin) density wave, which originates from alternating chains 
of Fe^"*" ions on octahedral Bl sites and Fe^^ ions on B2 sites (see Table I and Fig. 2 in 



Ref. 



471 ). A secondary [00|]c modulation in the phase of CO, which is formed by the chain 



of alternately "occupied" Fe^^ ions on the B4 sites and "empty" Fe'^"'" ions on B3 sites, was 
found. This is consistent with a [001] nesting vector instability at the Fermi surface in the 
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FeB minority electron states which has been recently revealed by the LSDA calculations for 
the cubic phase.— The calculated CO scheme coincides with the class I CO model proposed 
by Wright et al.—^ All the tetrahedra formed by Fe^ cations have either a 3:1 or 1:3 ratio 
of Fe^"*" and Fe^+ ions. Thus, the LSDA+f/ calculations confirm that the Anderson criterion 
is not satisfied in the LT phase. However, it should be pointed out that the Anderson 
criterion was introduced under the assumption of equal interatomic distances within each 
tetrahedron, whereas in the distorted LT structure the iron-iron distances vary from 2.86 to 
3.05 A. The same CO pattern has been recently confirmed by other LDA+f/ calculations.*^- 
An analysis of the 3d minority occupation matrices of Fee cations confirms very effective 
charge disproportion within the Fe^ minority spin subshell. In particular, one of the 
t2gi states of Fe^"^ and Fe^^ cations is almost completely filled with the occupation n ^ 0.8. 
On the other hand, the other two t2gi orbitals of the Fe^"^ cations have significantly smaller 
population of about 0.04. The occupation numbers of t2gi orbitals for Fe^"^ and Fe^ij cations 
do not exceed 0.1-0.17, which gives a value of about 0.7 for the largest difference of the 
populations of Fe^"*" and Fe^"*" t2gi states. The occupation numbers of the minority spin 
FeB 3d orbitals and the net occupations of the t2gi and Cgi states are given in the last two 
columns of Table [B 

The change of the t2gi occupations caused by the charge ordering is very effectively 
screened by the rearrangement of the other Fe electrons. A significant contribution to the 
screening charge is provided by Fe^ Cg states. Although the bands originating from these 
states are located well above the energy gap, the minority spin Cg orbitals form relatively 
strong a bonds with 2p states of the oxygen octahedron and, as a result, give an appreciable 
contribution to the occupied part of the valence band. The energy of Fe^"*" e^j states is 
lower than the energy of corresponding Fe^"*" states and the former give a significantly larger 
contribution to the part of the valence band formed mainly by O 2p states. Because of the 
stronger covalency of the Fe^"*" e^- O p bonds the net occupation of Fe^"^ Cgi states becomes 
~0.25e larger (see the last column of Table [11). The resulting 3d charge difference (0.23) 
and disproportionation of the total electron charges inside the atomic spheres of Fe^"*" and 
Fe^"^ ions (0.24) are in reasonably good agreement with the value of 0.2 estimated from a 
BVS analysis of the P2/c structure.— The above-mentioned screening of the changes in 
the Fe^ t2g minority occupations reduces the energy loss due to the development of charge 
order incompatible with the Anderson criterion in the LT phase of Fe304. 
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Hence, due to the strong screening effects, the order parameter defined as the difference of 
the net 3d charges of Fe^ cations does not provide conclusive evidence for CO. This explains 
why the BVS analysis does not give a convincing proof of CO existence. Apparently, a well- 
defined order parameter is the difference of the occupations of the minority spin states 
for Fe^^ and Fe^"*" cations which amounts to 70% of the ideal ionic CO model and clearly 
indicates the existence of a charge-ordered ground state below the Verwey transition. 

The LSDA+f/ calculations were also performed for the assumption of Verwey charge 
order in the P2/c structure. However, instead of the assumed Verwey CO the same self- 
consistent solution as the one described above was found. Therefore, the Verwey CO model 
is unstable in the distorted P2/c structure. It is well known that with increasing U value 
localization is effectively increased. Remarkably, even for a U value increased up to 7-8 eV 
no Verwey-like CO pattern was found self-consistently in the distorted P2/c structure. On 
the contrary, the LSDA+t/ calculations performed for an undistorted P2/c supercell of the 
Fd3m structure result in an insulating CO solution which is compatible with the Verwey 
CO model. Altogether this implies that the Verwey CO model is unstable under a structure 
distortion from the high-symmetry cubic into the low-symmetry P2/c phase. 

Also we performed LSDA+f/ calculations with the same U and J parameters (5 and 1 eV, 
respectively) for the assumption of one of the 32 class II CO models within Cc supercell of 
P2 / c, which is shown in Fig. 2 in Ref . |45|. But we found that this kind of CO is unstable and 
the self-consistent solution coincides with the one found previously for the P2/c structure. 

Comparing the LSDA+?7 results for the undistorted and distorted P2/c unit cells we can 
conclude that the charge-ordering pattern of Fe^"^ and Fe^"^ cations in the LT phase of Fe304, 



derived from the BVS analysis in Ref. |45| and confirmed by our study, is mainly forced by the 
local distortions of the crystal structure. Obviously we did not manage to study all possible 
charge- ordering scenarios within P2/c or Cc supercell of P2/c. But our results consistently 
indicate the importance of the small amplitude of atomic displacements (almost of 0.07A) 
recently resolved by x-ray and neutron powder diffraction.— The additional displacements 
leading to the Cc supercell were estimated to be of ~0.01 A but have not been fully resolved 
so far. They also may be important for full understanding of the CO in Fe304. In particular, 
in the P2/c subcell the true atomic positions are averaged over the corresponding number 
of subsites in the Cc cell. Therefore, the actual arrangement of the locally Fe^Og octahedra 
in the true Cc structure can be more complex, probably resulting in a more complicated 
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TABLE I: Sd orbital contribution to the formation of Fe^ minority spin states with occupancy n 
evaluated by diagonalization of the occupation matrix. Although one of the t2gi states of i?l and 
B4: sites is almost occupied with n > 0.7 e the t2g minority spin occupancies of B2 and B3 Fe^+ 
cations are less than 0.1 e. The occupied t2gi states of Bl and B4 Fe^"*" cations are predominantly 
of dxz^dyz and d^2_y2 character, respectively. The sum of t2g {eg) occupations is given in the last 
column. 
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charge and/or orbital order for the LT structure. The present calculations indicate that the 
competition of the "elastic" and electrostatic energy contributions in the total energy appears 
to be responsible for the CO, which is realized in the LT structure of Fe304.— Because of 
this, the Verwey CO model, which has the lowest electrostatic but significant "elastic" energy 
contribution in the total energy, becomes less favorable than other arrangements. 

C. Orbital ordering 

The self-consistent solution obtained by the LSDA+f/ is not only charge but also or- 
bitally ordered. It is clearly seen from Table [H which presents the contribution of 3d cubic 
harmonics to the formation of Fe^ minority spin states with an occupancy n (next to last 
column in Table [T]) evaluated by diagonalization of the corresponding occupation matrix 
self-consistently obtained by the LSDA+f/. 

As shown in the table the most occupied Fe^^ 3d minority orbitals are centered on the 
Bla, Bib, and B4 iron sites and have dxz — dyz, dxz + dyz, and d^^i-y^ character, respectively. 
Remarkably, the occupied t2gj, orbitals of Fe^ cations are almost orthogonal to each other, 
i.e., their relative orientation corresponds to an anti-ferro-orbital order. Since all Fe^ cations 
are ferromagnetically coupled the obtained orbital order conforms with the anti-ferro-orbital 
ferromagnetic state, which is the ground state of the degenerate Hubbard model according 
to the Kugel-Khomskii theory.^^'^^ This orbital order is consistent with the corresponding 
distortions of FeOe octahedra. In particular, using simple considerations which take into 
account only the change of the Fe-0 bond lengths and neglect the bending of the bonds, 
it was previously concluded that the calculated orbital order is mainly determined by the 
distortions of oxygen octahedra surrounding Fe^ sites.— 

Also this simple analysis shows a remarkable difference between Fe^"*]^ and Fe^"^ cations; 
namely, the average Fe^ia-O distance in the plane of occupied d^z — dyz orbital is 2.087 A, 
whereas in the planes of two other t2g orbitals they are only 2.063 and 2.067 A. This difference 
between the average cation-anion distance in the planes of occupied and unoccupied orbitals 
is remarkably larger for ^^Biafi (more than 0.02 A), although for FeB4 they are 2.074 and 
2.067 A for occupied d^^^yi and unoccupied d^z ± dyz orbitals, respectively, which gives a 
difference of only 0.007 A. This small difference can be changed by applying a uniaxial stress 
to the P2/c unit cell resulting in modification of the electronic properties. In particular, a 
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few percent of magnitude elongation of the P2/c unit cell along the c axis with simultaneous 
(in order to preserve the same unit cell volume) compression in the ab plane, gives rise to 
orbital-order crossover on the FeB4 site from a (i^.2_y2 to a d^z occupied orbital. At the 
same time the charge order and occupied orbitals on the Fe^^^^^ sites remain the same. 
The pressure-induced spatial reorientation of the occupied Fcba orbital was proved by 
the LSDA+f/ calculations for the strained P2/c unit cell. Note, however, that these rough 
estimations do not take into account the elastic anisotropy in Fe304. Moreover, the analysis 
was performed for the "averaged" P2/c structure. However, they provide insight into the 
orbital-ordering phenomena behind the Verwey transition in magnetite as well as the problem 
of an external parameter-controlled electron state (for example, orbital ordering) in solids.— 

D. Magnetic moments 

The strong variation of the occupancies of the minority spin Fe^ states leads to a 
pronounced modulation of the spin magnetic moments on the B sublattice. While the total 
moment per formula unit remains at 4/iB the magnetic moments of the Fe^^ Bl (3.50 ;Ub) 
and B4: (3.48 /xb) cations become appreciably smaller than FeB2 (3.94 /^b) and Fe^a (3.81 /^b) 
moments. The [001]c charge and spin modulation on the B sublattice is accompanied by 
formation of a weak spin modulation on the oxygen ions caused by different strengths of 
the hybridization of O 2p states with the minority spin 3d states of Fe^"^ and Fe^"*" ions. In 
particular, the oxygen magnetic moment reaches its maximal value of ^0.1 for 03 and 
04 ions, which lie in the plane of Fe^"^ cations. It substantially decreases for other oxygen 
ions and approaches minimum for 01 and 02 lying in the plane of Fe^'*[ cations (~0.04/iB). 

Recently, an anomalously large value of the Fe^ orbital magnetic moment reaching 0.33 fi-Q 
has been deduced by applying sum rules to experimental ^2,3 x-ray magnetic circular dichro- 
ism spectra of Fe304.'il In addition, the unquenched Fe^ orbital moment was also reported 
to be confirmed by the LDA+U calculations. Later, however, this experimental finding was 
questioned by Goering et al— The average orbital moments between —0.001 /xb and 0.06 /xb 
were found from x-ray magnetic circular dichroism (XMCD) sum rules depending on the 
integration range. From our spin-polarized relativistic LSDA+f/ calculations for the LT 
structure we obtained the orbital moments of 0.19/iB and 0.014 /^b for Fe^i and FeB2 ions. 
Somewhat larger values of 0.039 /ib and 0.22 /ib were calculated for Fbbs and Fsba cations. 



17 



respectively. Taking into account the negative Fe^ orbital moment of —0.021 /ie this gives 
the value of 0.07/iB for the average orbital moment. Thus, in agreement with the previous 



theoretical results of Ref. 



57 and XMCD sum rule data of Ref. 
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our calculations give the 



value of Fe orbital moment of ~0.07yUB which is much smaller than reported in Ref. 
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E. Exchange coupling constants 



B2(b) B3 




FIG. 5: (Color online) Sketch of the spatial arrangement of exchange interaction parameters 
between the octahedral Fe^ sites. Orbitals approximate the occupied 3d minority orbitals of Fe^ 
cations. Fe^ cations are shown by large (blue) spheres. Oxygen atoms are shown by small (green) 
spheres. Ferromagnetic couplings between Fe^ cations are shown by the thick (red) lines, whereas 
antiferromagnetic exchanges depicted by the thin (blue) lines. 



To proceed further we performed calculations of the exchange interaction parameters 
Jij via the variation of the ground state energy with respect to the magnetic moment 
rotation angle.— The exchange coupling parameter Jij represents the effective pair ex- 
change interaction between the ith and jth Fe atoms with effective Heisenberg Hamiltonian 
H = — J2i>j Jij^i ■ Sj. Here, Sj and Sj are the spins at site i and j (5/2 and 2 for Fe^+ and 
Fe^^ cations, respectively). Positive (negative) values of Jij correspond to the ferromagnetic 
(antiferromagnetic) coupling between sites. As shown in Table |TT] the exchange couplings 
between A and B iron sublattices are rather large, of about —70 K, and antiferromagnetic. 
The Fe^-Fe^ interactions are weakly antiferromagnetic with the maximal absolute value of 
9.3 K (not shown in Table UTl) . The exchange couplings between the Fe^ sites (| Jb_b| < 27.5 
K) are substantially smaller than Fe^-Fes ones and almost all of them are ferromagnetic (see 
Fig. [5]). Weak antiferromagnetic couplings with \ Jbb\ ^ 11-6 K are also obtained [mainly 
between the sites with the same 2+ or 3+ valence state, shown by the thin (blue) lines in 
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TABLE II: Exchange couplings Jij (all with \Jij\ > 10 K) are presented. The values are 
given in kelvin. The spatial representation of the Fes-Fes exchange couplings is schematically 
shown in Fig. [SJ Jij were calculated between the sublattices formed by the translations of 
the following Fe sites: Fe^i (1/4,0.0034,0.06366), FeA2 (1/4,-0.4938,0.18867), Fe^ia (0,1/2,0), 
FeBu (1/2,1/2,0), FeB2a (0,0.0096,1/4), Fessfe (1/2,0.0096,1/4), Fees (-1/4,0.2659,0.1198), 
Fe'53 (1/4,-0.2659,-0.1198), Fe'^3 (1/4,0.2659,0.3801), FeB4 (1/4,0.2479,-0.1234), Fe'54 
(-1/4,-0.2479,0.1234), and Fe'^4 (1/4,-0.2479,0.3765). 
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FIG. 6: (Color online) Comparison of the total DOS obtained from LSDA (dashed line) and 



LSDA+[/ calcu 
taken from Ref. 



ations with U =5 eV and J=l eV (solid line) to the experimental XPS spectrum 



361 (triangles). A rigid shift of —0.3 eV is applied to the calculated DOS curves. 



Fig. [5] . The spatial representation of these exchange couplings is presented in Fig. [51 Other 
couplings that are not shown in Table [III are weaker than 10 K. 

Experimental estimation of the exchange couplings in Fe304 was first performed by 
Neel on the basis of the two-sublattice coUinear model.— From analysis of the tempera- 
ture behaviour of the saturation magnetization and paramagnetic susceptibility he obtained 
Jaa = —17.7, Jab = —23.4, and Jbb = 0.5 K, where A and B refer to the tetrahedral 
and octahedral Fe sites, respectively. We find these values qualitatively in accord with our 
results presented in Table [Til namely, as in Neel's model, the calculations result in strong 
antiferromagnetic coupling between the A and B sublattices; Jaa couplings (not shown in 
Table [ITl) are considerably smaller than Jab] the exchange couplings in the B sublattice are 
weak and almost all of them are ferromagnetic. On the other hand, the small antiferro- 
magnetic Fe^'^^-Fe^'^^ exchange interaction (see Table [Til) is in exact agreement with recent 
estimations using the two-sublattice model. Three-sublattice model calculations give an 
overall similar result, except, however, the Fe^'^-Fe^'*' exchange coupling, which seems to be 
overestimated.— 
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VI. SPECTRAL PROPERTIES 



A. XPS spectra 

In Fig. [H] the LSDA and LSDA+f/ total DOSs are compared to the soft x-ray photoe- 
mission (XPS) spectrum from Ref. [SGj. No attempts to account for different cross sections 
of Fe 3d and O 2p states have been made. To account for the experimental resolution and 
finite lifetime broadening the theoretical spectrum was convoluted with a Gaussian with 
a half-width of 0.4 eV. The calculated DOS curves had to be shifted by 0.3 eV to lower 
energies in order to align the experimental peak at ~ —0.7 eV with the corresponding peak 
in the LSDA+f/ DOS. The necessity of the rigid shift can be attributed to the problem of 
experimental determination of the Fermi level position for insulating samples. A similar shift 



of 0.2 eV was used in Ref. |76| when comparing photoelectron spectra to LSDA calculations. 
The experimental spectrum was measured at 100 K using Fe 2p-3d resonance photoemission 
with the photon energy hu = 707.6 eV. Because of the the relatively large escape depth of 
photoelectrons at this photon energy it is expected that the spectrum mainly probes bulk 
electronic states of Fe304. 

Figure [6] shows that the shift of the position of the occupied majority spin Fes and 
minority spin Fe^ states to lower energies in the LSDA+[/ calculation greatly improves the 
agreement with the experiment. A feature at ~ —0.7 eV, which transforms to a well defined 
peak in the LSDA+f/ calculations, can be identified as — * d^ transitions at Fe^^ ions. This 
peak is completely spin polarized in accordance with the finding of a recent spin-resolved 
photoemission study. ''^ From a comparison with the /-projected DOS shown in Fig. [Hone 
can associate an experimental peak at 2 eV with transitions from the majority spin Fe^"^ Cg 
states which give the predominant contribution in this energy range. In the calculated DOS, 
however, this feature is masked by a peak at ~ —3 eV and appears only as a shoulder. A 
peak at 4 eV in the experimental spectrum originates to a great extent from the transitions 
from the majority spin Fe^"*" t2g strongly hybridized with O 2p states. Note, however, that 3d 
states of other Fe ions as well as O 2p states also contribute to the total DOS in the energy 
range from —7 to —4 eV. Finally, a DOS peak at —7.8 eV formed mainly by the majority 
spin Fe^"*" and minority spin Fe^ 3d states accounts for a weak feature found at —8 eV in 
the experimental spectrum. This feature seems to be masked by a high-energy satellite but 
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can be clearly seen in Fe 3p-3d resonance photoemission {hu = 56 eV).— 



B. Optical and magneto-optical spectra 

Motivated by the fact that previous calculations were performed for the undistorted 
cubic structure assuming the Verwey CO pattern,'^ in this section we present the optical 
and magneto-optical spectra calculated for the P2/c structural model of the LT phase of 
Fe304. Figure [7| shows optical conductivity (top) and reflectivity (bottom) calculated using 
the LSDA and LSDA+f/ with U=5 eV and J=l eV. Symmetry lowering from the fee to P2/c 
structure causes some degree of optical anisotropy in the low-frequency range. However, it 
is rather weak and here we present the spectra averaged over different photon polarizations. 
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The results are compared to the experimental optical spectra taken from Ref. 
reveal strong temperature dependence only in the narrow energy range of 0-1 eV. 

The LT structure distortions have a rather weak effect on the LSDA spectra which re- 
semble closely the spectra calculated for the cubic phase.— This finding agrees well with the 
fact that the LSDA band structure presented above is not strongly affected by the crystal 
structure distortions. It is worth recalling that the LSDA calculations result in a metallic 
solution and therefore give the wrong asymptotic behavior for the optical reflectivity at 
uj —>■ 0. Moreover, a peak at ~6.5 eV in the calculated conductivity spectrum is shifted to 
higher energies with respect to the experimental peak centered at ~5 eV. Both conductivity 
and reflectivity show peaks at 2 eV which are much higher than in the experimental spectra. 
Hence, the LSDA theory fails to correctly reproduce the low-energy spectral properties of 
Fe304, whereas it describes reasonably well the experimental spectra above 10 eV. 



As in Ref. |57| better agreement between the calculated and experimental spectra was ob- 
tained using the LSDA+f/ method. Since the self-consistent LSDA+t/ solution is insulating 
the reflectivity at zero frequency is no longer equal to 1 but approaches the value of 0.6 
which, however, is somewhat higher than in the experiment. Figure [8] shows the expanded 
view of the theoretical and experimental^ optical spectra in a smaller frequency range below 
8 eV. Charge and orbital ordering within the minority spin Fe^ states changes drasti- 
cally the spectrum below 3 eV. Because of the shift of the unoccupied minority spin Fe^"^ 
and majority spin Fe^ d states to higher energies as compared to their LSDA positions the 
spectral weight is transferred from 2 eV to higher frequencies. This improves the agreement 
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FIG. 7: (Color online) Reflectivity (lower panel) and optical conductivity (upper panel) spectra 
calculated using the LSDA (dashed line) and LSDA+[/ with U=5 eV and J=l eV (solid line) for 
the P2/c model of the LT phase of Fe304. Room temperature experimental spectra from Ref. 
are shown by triangles. 
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FIG. 8: (Color online) Comparison of LSDA (dashed line) and LSDA+[/ with U=5 eV and J=l eV 
(solid line) optical conductivity spectra to the experimental spectra measured at room temperature 
(up triangles) and at 10 K (down triangles).— Also plotted are the contributions of the interband 
transitions from occupied minority spin Fe^"*" t2g to Fe^"*" t2g (short-dashed line) and Fe^"*" eg (dash- 
dotted line) to the LSDA-I-C/ spectrum. 
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with the experimental optical conductivity in the range 3.5-6 eV. However, instead of the 
two-peak structure at 0.5 and 2 eV the LSDA+[/ calculations reveal a maximum at ~1 eV 
which can be associated with the 0.5 eV peak in the experimental spectrum. The theoretical 
peak is, however, more symmetric, shifted to higher frequencies, and its magnitude is lower 
than in the experiment. An analysis of different interband contributions to the calculated 
spectra shows that the 1 eV peak is predominantly formed by the transitions from occupied 
Fe|"^ t2gi into empty Fe^"*" t2gi bands. The maximum of the Fe|"^ t2gj, to Fe^"^ e^j interband 
conductivity coincides with the second experimental peak at ~2 eV. However, its theoret- 
ical intensity is too low and these transitions do not manifest themselves as a peak in the 
spectrum. 

An attempt to make the agreement between the theoretical and experimental spectra 
better by decreasing the value of U has not given any substantial improvements. In particular 
for f/=4.5 eV (not shown here) the 1 eV peak shifts to lower energies only by 0.2 eV. 
Moreover, its shape is still more symmetric than the shape of the corresponding peak in the 
experimental spectra at 10 K,^^ whereas the theoretical intensity of the second peak remains 
too low. 

A possible reason for the abovementioned discrepancies between the theoretical and ex- 
perimental spectra below 2 eV can be neglecting dynamical correlations in the frame of the 
static LSDA+[/ method. Moreover, it is worth recalling that the true LT crystal structure 
has not been fully resolved yet and the calculations were performed for the "averaged" P2/c 
structural model for the LT phase of Fe304. As described in Sec. IV Bl the self-consistently 
obtained CO corresponds to the class 1 CO model.— If one of the class II CO patterns is 
realized in the true Cc structure, as has been suggested by a recent resonant x-ray diffrac- 
tion study,— then more complicated charge and orbital order may result in narrowing of Fe^ 
d-dehved bands and a shift of the peak in optical absorption to lower energies. Therefore, 
further investigations of the LT crystal structure are strongly demanded. 

In Fig. [9] magneto-optical polar Kerr rotation and ellipticity spectra calculated for the 
P2/c structural model are compared to the room temperature spectra of the annealed syn- 
thetic sample of Fe304.™ Calculations were performed assuming that the magnetization and 
light propagation directions are parallel to the c axis of the P/2c cell. Similar to the optical 
conductivity spectra the use of the LSDA+f/ approach improves the agreement with the ex- 
periment above 2 eV. The minimum at 1.3 eV in the theoretical Kerr rotation spectrum is. 
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FIG. 9: (Color online) Polar Kerr rotation (lower panel) and ellipticity (upper panel) calculated 
for the P2/c model of the LT phase of Fe304 using LSDA (dashed line) and LSDA+C/ with U=5 
eV and J=l eV (solid line). Experimental spectra measured at room temperature for an annealed 



synthetic sample of Fe304 in Ref. 



79 are plotted by triangles. 



however, less pronounced and shifted to higher frequencies as compared to the experimental 
one (measured at room temperature). Comparison of the magneto-optical spectra for the 
LT phase to the spectra calculated assuming the Verwey CO in the fee unit cell^ shows 
that the theoretical spectra are rather sensitive to the structural and CO model used in 
the calculations. Its experimental measurements for the LT phase of Fe304 and comparison 
with the theoretical prediction presented in Fig. M are of strong importance. 



C. O K-edge XAS 

In the process of x-ray absorption at the O K edge a core O Is electron is excited into 
unoccupied states in the conduction band. In the dipole approximation, which was used 
in the present calculations, the final states for the transition are of p symmetry and the 
spectrum probes the density of unoccupied O 2p states. O i^'-edge x-ray absorption spectra 
calculated using the LSDA and LSDA+f/ (with U=5 eV and J=l eV) and averaged over 
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FIG. 10: (Color online) LSDA (dashed line) and LSDA+J/ with U=5 eV and J=l eV (sohd line) 
O K-edge x-ray absorption spectra calculated for the P2/c model of the LT phase of Fe.s04. The 
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8C 



sum of the spectra from inequivalent O sites is shown. Experimental spectra from Refs. 
are plotted by down and up triangles, respectively. 

eight inequivalent O sites in the P2 / c unit cell are presented in Fig. [TOl Experimental spec- 



tra taken from Refs. 
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and [801 are shown in Fig. by down and up triangles, respectively. 
In order to account for the effect of the final life time of the O Is core hole and experimental 
resolution the spectra were broadened by convolution with Lorentzian and Gaussian func- 
tions of the width of 0.4 and of 0.5 eV, respectively. Then the calculated spectra were aligned 
and normalized to the main absorption peak at 541 eV in the experimental spectrum. 

As mentioned in Sec. |IV]the lattice distortions themselves have only a weak effect on the 
LSDA band structure of Fe304 and the LSDA spectra calculated for inequivalent O sites in 
the P2/c cell resemble closely the O i^'-edge spectrum for the undistorted Fd3m structure. 
Both the LSDA and LSDA+[/ spectra reproduce well the shape of the 541 eV peak and the 
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that these 



position of a peak at 547.5 eV. Our calculations confirm the conclusion of Ref. 
peaks originate from transitions to bands formed mainly by delocalized Fe 4s and Ap states. 

The so-called preedge peak at 531 eV is formed by transitions to O 2p states hybridized 
with FeA and Fe^ 3d states and, thus, is more sensitive to the changes of the electronic 
structure caused by charge and orbital ordering effects. The preedge peak in the LSDA 
spectrum is wider than the experimental one, more intense, and shifted to lower energies. 
The experimental peak exhibits a second feature that appears as a shoulder at the high- 
energy slope of the peak. The LSDA spectrum also has a two-feature structure but their 
intensities are reversed. The two-feature structure of the preedge peak was interpreted as 
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coming from the hybridization of O 2p states with Fe 3d states spht by the hgand field.— 
Indeed, the comparison of the O K-edge spectrum calculated for the undistorted fee structure 
to the corresponding DOS shows that the low-energy shoulder is formed by transitions to 
O 2p states hybridized with majority spin FeA 3d states of e symmetry and minority spin 
FeB states, whereas the main peak reflects the O 2p hybridization with Fe^ ^ind Fe^ 
Cgi states. In the latter case the O 2p - Fe 3d hybridization is stronger and this causes the 
higher intensity of the corresponding feature in the LSDA spectrum. It should be noted that 
calculations of the O K-edge spectrum based on the multiple-scattering approach reproduced 
the relative intensities of the two features of the preedge peak.— 

O -ft"-edge spectra calculated for inequivalent O sites in the P2/c unit cell using the 
LSDA+f/ approach show much stronger variation of the shape of the preedge peak as com- 
pared to the LSDA calculations. The strongest deviation of the peak shape is found for 
those O ions that have only Fe^"*" cations among their nearest neighbors. In spite of that, 
the preedge peak in the averaged spectrum becomes narrower than the LSDA one and shifts 
to higher energies. An additional feature does appear at the high-energy slope of the peak 
but its intensity is much weaker than in the experiment. 

We can conclude that LSDA+f/ calculations reproduce better the energy position of the 
preedge peak in the O K-edge spectrum but the shape of the peak still differs from the 
experimentally observed one. Calculations performed for other CO models possible in the 
true Cc structure might in principle improve the agreement with the experiment. We have 
to note, however, that recent measurements of the O i^-edge spectrum performed in a wide 
temperature range below and above the Verwey transition show very small variation of the 
shape of the preedge peak.— This implies that the O K-edge spectrum should not be very 
sensitive to the details of a particular CO pattern. Many-body effects that are beyond the 
LSDA+f/ may be responsible for the disagreement between the theory and the experiment. 

VII. SUMMARY AND CONCLUSIONS 

In summary, in the present LSDA+f/ study of the P2 / c model of the LT phase of Fe304 we 
found a charge- and orbitally ordered insulator with an energy gap of 0.18 eV. The obtained 
charge-ordered ground state is described by a dominant [001]c charge density wave with 
a minor [00|]c modulation on the Fe^ sublattice. A weak [001]c spin/charge modulation 
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on the oxygen ions was also obtained. The CO coincides with the earher proposed class 
I CO^ii^ and confirms violation of the Anderson criterion.^- While the screening of the 
charge disproportion is so effective that the total 3d charge disproportion is rather small 
(0.23), the charge order is well pronounced with an order parameter defined as the difference 
of t2gi occupancies of 2+ and 3+ Fe^ cations (0.7). This agrees well with the result of 
BVS analysis for a monoclinic structure (0.2). The orbital order is in agreement with 
the Kugel-Khomskii theory®^ and corresponds to the local distortions of oxygen octahedra 
surrounding Fe^ sites. The average Fe orbital moment of ~ 0.07 agrees well with the 
recent experimental findings.— 

Calculations of the effective exchange coupling constants between Fe spin magnetic mo- 
ments show that the dominating interaction is an antiferromagnetic coupling between Fe^i 
and Fe^ moments. The coupling between Fe^"*" and Fe^"^ moments is found to be weaker 
and ferromagnetic. 

The relevance of the charge-ordered LSDA+f/ solution has been verified by performing 
calculations of optical, magneto-optical polar Kerr rotation, and O i^-edge x-ray absorption 
spectra and comparing them to available experimental data. The LSDA+f/ spectra show 
much better agreement with the experimental ones as compared to the spectra calculated 
for the charge-uniform LSDA solution. Unfortunately, the remaining discrepancies between 
the theory and the experiment in the low-frequency part of the optical and MO spectra and 
in the shape of the preedge peak in x-ray absorption do not allow us to discriminate among 
different CO models in the true Cc structure. 
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